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Using high-speed photography coupled with optical interference, we experimentally study 
the air entrapment during a liquid drop impacting a solid substrate. We observe the 
formation of a compressed air film before the liquid touches the substrate, with internal 
pressure considerably higher than the atmospheric value. The degree of compression 
highly depends on the impact velocity, as explained by balancing the liquid deceleration 
with the large pressure of compressed air. After contact, the air film expands vertically 
at the edge, reducing its pressure within a few tens of microseconds and producing a 
thick rim on the perimeter. This thick-rimmed air film subsequently contracts into an air 
bubble, governed by the complex interaction between surface tension, inertia and viscous 
drag. Such a process is universally observed for impacts above a few centimeters high. 



1. Introduction 

The impacts of liquid drops onto solid substrates are ubiquitous and appear in a variety 
of applications, such as spray coating and ice accumulation on aircraft ( Yarinll2006f ). While 
the interaction between liquid and solid has been extensively studied, the important role 
of air is discovered only recently, as shown b y the surprising finding t hat air pressure 
strongly influences the liquid splash outcomes ( Xu et a/1 120051 : IXul 12007 ). Therefore, un- 



derstanding the behavior of air during liquid-solid impacts will bring new advancement 
to this fundamental phenomenon, and may benefit practical processes such as splash con- 
trol and surface coating. Experimentally, entrapment of air during impacts is commonly 
observed: when the impact speed is high (~ lm/s), a thin air film is t rapped at the very 
beginning, which subsequently contracts into one or two air bubbles ([Thoroddsen et all 
20031 [2QQ5[ ). At low impact speed (~ O.lra/s), however, the air under the drop remains 
connected with outside for the majo rity of time (i.e., not entirely trappe d), until it gets 
enclo sed by the moving contact line (jde Ruiter et al\\2012\) . Recent study (jKolinski et al. 



20121 ) further reveals a short-lived nanometer thick air film right before the formation 
of contact line. The air entrapment depends on the impact velocity, the liquid property 
([Thoroddsen et aZ.I I2QQ31 120051: Ide Ruiter et al\\20l£ iKolinski oZll2Q12h and the sur- 
rounding air pressure ([Driscoll fc N agel 2oTl|). Theoretically, simulations have explored 
air entrapme nt with both compressible dMandre et al. 20091 Mani et ~al\ l2Q10h and in- 
compressible ( Smith et a/.ll2003[ iHicks fc Purvisll2010l ) models. Marked deformation of 
drop surface before it touches the solid is predicted: the drop surface is deformed upward 
at the impact center, m aking the first c ontac t awa y from the c e nter o n a ring-like area. 
Compressible model by iMandre et all ( 20091 ) and iMani et a/1 ( 2010l ) further indicates 
significant compression in the trapped air, resulting in a pressure considerably higher 
than the atmospheric value. Despite these important predictions, however, experimental 
measurement on the exact condition of trapped air is still missing, especially near the 
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liquid p(kg/m 3 ) a(mN/m) /x(mPa • s) ^(degree) R(mm) 

H 2 1000 50 ±5 1.00 ±0.01 65 ±5 2.0 ± 0.1 (1.5 ± 0.1) 

oil-1.04 816 17.4 1.04 ±0.01 1.5 ± 0.1 (1.2 ± 0.1) 

oil-9.30 930 20.1 9.30 ± 0.02 1.5 ±0.1 

Table 1. material properties and drop radius of different liquids. R is measured right before 
contact by fitting the local radius of curvature at the bottom of the drop. For the lowest release 
height, strong oscillation in the drop makes R different from the average radius in H2O and 
oil-1.04, as shown in the parenthesis. 



critical moment of impact. In particular, even the fundamental question whether the 
trapped air is compressed or not remains unclearly. To clarify these puzzles and better 
understand the impact phenomenon, experimental study on air entrapment close to the 
moment of impact is highly desirable. 



2. Experimental methods 

Using fast photography coupled with optical interference ( Driscoll fc Nagelll2Q 1 ll : Ide Ruiter et al. 



l2012f h we experimentally study the air entrapment during a liquid drop impacting a 
smooth substrate at relatively high speeds (0.7 — 3m/s). We observe the formation of a 
compressed air film before the liquid touches the substrate, with internal pressure con- 
siderably higher than the atmospheric value. The degree of compression highly depends 
on the impact velocity, which is explained by balancing the liquid deceleration with the 
large pressure of compressed air. After contact, the air film expands vertically at the 
edge, reducing its pressure within a few tens of microseconds and producing a thick rim 
on the perimeter. This thick-rimmed air film subsequently contracts into an air bubble, 
governed by the complex interaction between surface tension, inertia and viscous effects 
(see Fig. 5 (a) for the entire process). Such a process is universally observed for impacts 
above a few centimeters high. 

To independently study the effects of viscosity and surface tension, three different 
liquids are used: H2O, oil-1.04 and oil-9.30. The two oils are silicone oils with similar 
surface tensions a but different dynamic viscosities: \i = 1.04mPa • s and 9.30mPa • s. By 
contrast, H2O and oil-1.04 have similar viscosities but different surface tensions (see table 
[T]). Reproducible drops around millimeter in size (see table Q] for details) are released from 
rest at different heights, and impact a smooth and dry cover glass at various velocities 
(0.7 — 3m/s). The two silicon oils completely wet the glass substrate with zero static 
contact angle, and H2O has the static contact angle = 65° ± 5°. All experiments 
are performed at the atmospheric pressure, Po = lOlkPa. The impacts are viewed from 
below with an inverted microscope and recorded with a high-speed camera (Photron 
SA4) at recording speeds up to 150,000 frames per second. The illumination light is 
monochromatic, with the wavelength A = 546nm and coherence length a few microns. 
The short coherence length makes sure that there is no interference between the two 
sides of the substrate. The Newton's rings produced by the trapped air can quantitatively 
characterize the thickness profile of air. 
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Figure 1. (color online) The formation of trapped air film before liquid touches solid, (a) The 
interference patterns before contact by an oil-9.30 drop at impact velocity Vo = 0.74 ± 0.01 m/s. 
The time uncertainty is half of the frame interval, 4/xs. In the first frame (— 20/xs), Newton's 
rings appear as the drop-to-substrate distance becomes smaller than the coherence length of 
light. In the second frame ( — 13/xs), a dimple with horizontal scale several hundred microns but 
vertical scale a few microns is produced. The third frame (— 7/xs) reveals a few black spots of 
isolated contacts at the edge (specified by the white arrows). The fourth frame (0/is) shows the 
global liquid-solid contact on a thick black ring at the edge of the pattern. This moment of first 
global contact is also defined as the moment of t = throughout our study, (b) The thickness 
profile, H(r), constructed from the patterns in (a). Note the very different scales for x and y 
axes. The curves almost overlap with each other, indicating very little motion of the interface 
through different frames. Significant difference only appears at impact center (r = 0), where the 
interface moves upward instead of downward, with the speed 0.02±0.01m/s. The schematics in 
the inset illustrates a dimple with horizontal scale 2Lo and thickness Ho (drawing not on scale). 



3. Experimental results 

3.1. Formation of air film before contact (t < 0) 

We first clarify the entrapment behavior before liquid touches solid. As the drop ap- 
proaches the substrate, the air in the thin gap can not escape immediately and gets 
compressed to a high pressure (exact values measured later in Fig.4(b)). Such a high 
pressure locally pushes in the drop surface and creates a dimple around impact center, 
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as illustrated by the cartoon in Fig. 1(b) inset. The interference patterns from the dimple 
are shown in Fig. 1(a), for an oil-9.30 drop at the impact speed Vo = 0.74 ± O.Olm/s. In 
the first frame (— 20/xs), Newton's rings appear as the drop-to-substrate distance becomes 
smaller than the coherence length of light (a few microns). In the second frame (— 13/is), 
a quite flat dimple with the horizontal scale several hundred microns but the vertical 
scale a few microns is produced. Because of this thin and flat nature, the dimple can be 
exactly considered as an air film. The third frame (— 7/is) shows a very similar pattern 
as the second one, except with a few black spots of isolated contacts indicated by the 
white arrows. The fourth frame (0/is) reveals the global liquid-solid contact on a thick 
black ring at the edge of the pattern. This moment of first global contact is also defined 
as the moment of t = throughout our study. 

The thickness of the dimple, H, can be accurately derived from the pattern with the 
relationship AH = A/4 = 136.5nm between two neighboring rings. Using contact area as 
the zero-thickness reference point, we can quantitatively determine the entire thickness 
profile of the trapped air in the t = frame. We can further determine the profiles for 
the two previous frames (t = — 7/is and — 13/is), due to their highly similar patterns as 
the t = frame, which enables the tracking of every ring. The thickness as a function 
of distance from the center, H versus r, is plotted in Fig. 1(b): 13/is before the contact, 
an air film with the lateral size 2Lo = 500/mi but a thickness of only H = 3.5/mi forms, 
which maintains an almost identical profile through the next several frames. Significant 
variation only appears near the center (r = 0), where the interface moves upward instead 
of downward, with a small speed 0.02 ± O.Olm/s (calculated from the profiles at t = — 7/is 
and t = 0/is). Clearly, the motion of the drop is quite complex immediately before the 
contact: while the main body falls at the impact speed, Vo = 0.74m/s; the small volume 
above the trapped air (2Lo = 500/im) stays almost stationary, with the region at r = 
even moving oppositely in upward direction. 

3.2. Profile of air film upon contact (t = 0) 

To systematically study the air film properties, we characterize their profiles at t = 
for different velocities and liquids, as shown in Fig. 2 (a). The three panels correspond to 
the three different liquids and each panel contains data for several impact velocities. In 
general, air film becomes thinner with the increase of impact velocity, as illustrated by the 
measurements. For quantitative understanding, we plot the maximum thickness measured 
at r = 0, Ho, as a fun ction of Vp in Fig. 2(b ). Clearly Ho decreases with Vo, consistent with 
the previous theory ( Mandre et al. ( 2009f )). However, the exact dependence of Hq on Vo 



differs from the predictio n: we calculate the "c ompressible factor" defined theoretically 
as e = Pn/(R/i- 1 V^ 4 ) 1 / 3 dMaii(ire et al\ (j2QQ9h . fig is the dynamic viscosity of air), and 



obtain the range e ~ 0.1-10 for our experimental data. Thus our experiments are at 
the transition region from "incompressible regime" to "compressible regime" predicted 

2/3 

by the theory, and Ho should decrease with Vo faster than the power law of Ho ~ V 
However the log-log plots in Fig. 2(b) find powers much slower than -2/3. This discrepancy 
calls for further theoretical studies on the problem. 

We also measure the hori zontal radius , Lp, v ersus Vo for different liquids and compare 
it with the existing theory ([Hicks et all (2012)). In Fig. 2(c), we plot our m easurements 



as soli d symbols and the theoretical prediction, Lo = 3.8 • (4/x g /pVo) 1 ^ 3 R 2//3 ([Hicks et all 
(2012)), as the open symbols. We use the bottom radius of curvature right before contact 



shown in Table 1 for the calculation. Without any fitting parameter, a reasonable agree- 
ment is observed at high Vo while some deviation occurs for low Vo- This is surprising 
since the theory is based on an incompressible calculation, which should match our low 
Vo rather than the high Vo region. Further study is required to clarify this puzzle. 
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Figure 2. (color online) Characterizing the air film profile at t = for different liquids and 
impact velocities, (a) Air film profiles measured at the moment of contact (t = 0) with New- 
ton's rings method. The three panels show the results for H2O, oil- 1.04 and oil-9.30 respectively. 
Different colored curves are measured at different Vo- (b) The maximum film thickness at the 
center, Ho, v.s. the impact velocity, Vo- The solid lines are the best power- law fits within the 
limited data range, (c) The horizontal radius, Lo, versus impact velocity, Vo, in different liquids. 
The solid symbols are experimental measurements, while the open symbols are theoretical pre- 
dictions from an incompressible model, L = 3.8 • (4/i g /pVo) 1/3 R 2/3 (Hicks TFall (|2012h ). Good 
agreement is observed at high Vo. 



3.3. Evolution of air film after contact (t ^ 0) 

The dimple continues to evolve after the liquid-solid contact. The interference patterns for 
t ^ are shown in Fig. 3 (a), for an impact with the same condition as in Fig. 1(a). Frame 
one (0/is) shows the Newton's rings from the initial thin film, with the global contact 
on a thick black ring at the edge. During the subsequent contraction (28 — 139/is), the 
Newton's rings in the middle remain largely unchanged, while a grey region without 
rings grows thicker and thicker at the edge. This grey region smoothly evolves into an air 
bubble at the end of the process (319/iS and 2486/is). Such a smooth evolution implies 
that the thickness of grey region at t = 319/is must be close to the bubble diameter 
(d = 76/im) and substantially exceed the coherence length of light, which explains the 
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Figure 3. (color online) The evolution of trapped air after impact. The time uncertainty is 
half of the frame interval, 7 lis. (a) The interference patterns for t ^ by an oil-9.30 drop at 
Vo = 0.74 ± O.Olm/s. Frame one (0/is) shows the Newton's rings from the trapped thin air film, 
with the thick black ring at the edge indicating global liquid- solid contact. During the subsequent 
contraction (28 — 139/xs), the Newton's rings in the middle remain largely unchanged, while a 
grey region without rings grows thicker and thicker at the edge. This grey region smoothly 
evolves into an air bubble with diameter d = 76 =b 3/rni (319/xs and 2486/xs). (b) The thickness 
profiles from the center towards the edge along one typical radial direction. From the patterns in 
(a), we derive the solid symbols using the Newton's rings, and estimate the grey region's profile 
with half circles as plotted by the dashed curves (Although different scales in x and y axes 
produce an non-circular impression). We emphasize that the dashed curves are from estimate 
and thus can not be taken as serious measurements. 



lack of Newton's rings. Similarly, having a thickness larger than the coherence length 
explains the absence of rings in the grey region for previous frames as well. To estimate 
the thickness there, we assume that the grey region has a half-circular vertical cross 
section, as demonstrated by the dashed curves in Fig. 3(b) (the non-circular impression 
comes from the different x and y scales). 

Combining the thickness measured from the Newton's rings and estimated at the grey 
region, we plot the entire thickness profiles for the t ^ frames in Fig. 3(b). Different 
colors indicate different time t, with the solid symbols from direct measurement and the 
dashed curves from the half-circular estimate. Clearly the air film does not contract with 
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Figure 4. (color online) Characterizing the initial pressure (i.e. at t = 0) for different liquids 
and impact velocities, (a) The volume of the initial air film, Vmm, and the final air bubble, 
Vbubble, v.s. the impact velocity Vo in different liquids. Apparently the difference between Vmm 
and Vbubble increases with Vo- (b) The pressure of air film at t = 0. The solid symbols come 
from volume measurements fEa !3 ip and the open symbols are calculated from the expression 
fEa !3.2p . They agree reasonably well. The fitting parameter, C, has the values: C = 0.33 (H2O), 
0.33 (oil-1.04) and 0.23 (oil-9.30). It does not change with surface tension but decreases with 
viscosity. 

a uniform thickness; instead the edge grows into a thick rim while the interior remains 
thin and fiat. 

3.4. The compression at t = 

Is the trapped air initially (t = 0) compressed or not? We address this fundamental 
question by accurately measuring and comparing the volumes at the initial (t = 0) and 
the final states. The volume at t = is directly obtained from the thickness profile H(r); 
and the volume of the final air bubble is accurately determined from its diameter. To 
make sure that the bubble is a perfect sphere instead of a cap intersected by the substrate, 
we keep track of the bubble until it leaves the surface and rises upward. Since the bubble 
is a perfect sphere after leaving the surface, we can unambiguously determine its volume 
from the diameter. The exact values of Vm m and Vbubble are shown in Fig. 4(a) for 
different velocities and liquids. The plot reveals a considerable volume increase from film 
to bubble, proving that the trapped air at t = is indeed compressed. 

To further estimate the pressure at t = 0, we need to clarify whether the compres- 
sion is isothermal or adiabatic. It depends on the time scale of air entrapment, r, com- 
pared with the time scale during which a thermal equilibrium can be reached, r f . In 
our experiment, r ~ Ho/Vo ~ while r f is determined by the rate of thermal 

conduction. Our system is a thin air film with the upper and lower boundaries (liq- 
uid drop and glass substrate respectively) at the room temperature. Suppose there 
is a temperature change AT in the middle (from gas compression or viscous heating 
or other reasons), the thermal gradient then becomes AT/(Hq/2) and the heat flux 
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per unit time is dQ/dt ~ 2 • k • A • AT/ (Ho/2), with k being the thermal conductiv- 
ity of air and A being the horizontal cross section area; the factor of 2 comes from 
the existence of two boundaries. This heat flux will bring the system back to equi- 
librium during the time scale: r' ~ c p • p a i r • v a i r • AT/(dQ/dt) ~ c p • p a i r • Ho/12k, with 
c p being the air's specific heat per unit mass and v a i r ~ ^A • H being the volume of 
trapped air. Plugging in the values for atmospheric pressure and room temperature, 
c p = 10 3 J/(kg-K),p air = 1.2kg/m 3 ,k = 0.026J/(s-m-K), we obtain: r' ~ 0.03/is <C 
r ~ lfis. Therefore, the trapped air reaches thermal equilibrium rather rapidly, due 
to its thin thickness and small heat capacity, and the compression in our experiment can 
be considered as isothermal. Subsequently, the compressed air film expands to an uncom- 
pressed air bubble, again in an isothermal manner since expansion is even slower than 
compression. Thus we can obtain the initial air film pressure, Pfiim, from the isothermal 
equation of state: 



-P, -p, ^bubble -p, Vbubble /o -i \ 

^film = -^bubble " 77 = r*0 • — (6.1) 

Vfiim Vfiim 

Here Pbubbie = Po = lOlkPa since the bubble is uncompressed (the curved surface of 
bubble adds a negligibly small Laplace pressure around lkPa to Pbubbie) - The exact values 
of Pfiim are shown as solid symbols in Fig. 4(b). Apparently, Pfii m increases dramatically 
with impact speed, varying from very close to Po to several Po- 

This velocity dependence can be qualitatively understood by a force balance estimate. 
Immediately before the contact, the compressed air in the dimple locally decelerates the 
liquid to an almost complete stop in vertical direction (see Fig. 1(b) main panel), on the 
length scale of the dimple's radius Lo (see Fig. 1(b) inset). The deceleration occurs dur- 
ing the time scale r ~ Ho/Vo, with the magnitude a ~ Vo/r ~ Vq/Ho. Since the affected 
liquid has the mass m ~ pLjj, we get the force: f = m • a ~ CpLjj Vq/Ho, with C being 
a pre-factor of order unity. This force naturally comes from the excess pressure of com- 
presse d air multiplying the area: f ~ (Pfi i m ~ P o) ■ Lq- Making the two force expressions 



pressed air multiplying tne area: 1 ~ (rf n m — r p ) • Lq. Making the two lorc( 
equal (iMandre et all (|2QQ9h : iMani et all <|2010h ) and solving for P fi i m yield: 



Pfiim = Po + C • pYl • ^ (3.2) 

Wo 

This expression has the same trend as the measurements: at small Vo, Pfiim approaches 
Po; while it rises rapidly with Vo due to the large geometric factor, Lo/Ho ~ 100 ^> 1. 
The values calculated from this expression are plotted as open symbols in Fig. 4(b), and 
agree reasonably well with the measurements (solid symbols), for different liquids at 
various impact speeds. The only fitting parameter, C, is indeed of order unity: C = 0.33 
(H2O), 0.33 (oil-1.04) and 0.23 (oil-9.30). Clearly C does not change with surface tension 
but decreases with viscosity, suggesting that under the same impact condition the air 

trapped in a more viscous drop is less compressed. 

Comparing with the existing compressible model (iMandre et all ( 20091 ): IMani et ~al. 



(|2Qinh ). our experimental range overlaps with their transition region from "incompress- 



ible regime" to "compressible regime" (i.e., 0.1 < e _1 < 10). Indeed the data in Fig. 4(b) 
verify such a predicted transition: at low Vo the pressure is close to Po and the compres- 
sion is small; while significant compression appears above Vo ~ 1.5m/s (corresponding 
to e ~ 1 in the theory). We also find that the initial impact pressure, Pfu m , exceeding 
the liquid inertia, pVg, by a large factor of Lo/Ho ~ 100. This comes from the dramatic 
deceleration of liquid under the thin gap geometry. 
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Figure 5. (color online) Evolution of trapped air during contraction, (a) Schematics illus- 
trating the vertical profile evolution during contraction. The air film turns into uncompressed 
thick-rimmed structure within a few tens of microseconds, and then contracts into an air bubble. 
(6) Measuring S(t) experimentally. The three curves come from the three liquids at the same 
impact velocity, Vo = 1.49 ± 0.02m/s. The open symbols are our measurements and the solid 
curves are high-order-polynomial fits. These functions are plugged into Eq.3.4 to calculate L(t). 
(c) The comparison of L(t) between the measurements and the calculation. The open symbols 
come from direct measurements while the curves are calculated from Eq.3.4. With the same 
fitting parameters, Ci = 0.63 and C2 = 1.54, we obtain excellent agreement for different liquids. 
The air film contracts the fastest in H2O, due to the large surface tension and low viscosity. It 
contracts slower in oil- 1.04, because of the reduced surface tension. The motion is the slowest 
in oil-9.30, due to the effects of both low surface tension and high viscosity. 




3.5. The contraction at t > 

From compressed air film to uncompressed air bubble, the detailed process requires fur- 
ther clarification. We measure this process from interference patterns and demonstrate 
the evolution in Fig. 5 (a). At the moment of contact (t = 0), the trapped air is thin and 
flat with a pressure higher than Pq. After the contact (t > 0), however, the air pressure 
drops sharply as the impact pressure decreases, resulting in a rapid volume expansion at 
the edge of film. This non-uniform expansion may come from the non-uniform pressure 
distribution after contact: the center of the film is right below the stagnation point of 
liquid and thus experiences higher impact pressure; while the edge of air film experiences 
less impact pressure and more easily to expand. This expansion produces a rim much 
thicker than the interior, as illustrated in Fig. 5 (a). We estimate the rim's volume by 
assuming a half-circular profile and deduce the pressure from volume estimate. We find 
that the pressure rapidly drops to uncompressed value within a few tens of microseconds. 
The uncompressed thick-rimmed structure subsequently contracts into an air bubble. 

Therefore, the contraction process is mostly under uncompressed condition except at 
the very beginni ng, consistent with the incomp ressible assumption in previous contraction 
measurements ([Thoroddsen et q/J I20Q3L l2005h . However, the previous study assumes a 
uniform film thickness, while we believe that the thick rim plays an important role. We 
take the characteristic size of the rim, S (see Fig.5(a)), as the dominant length scale during 
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contraction: the surface tension provides the driving stress, a/S, which is balanced by the 
liquid inertial, pv 2 , and the viscous stress, fiv/5. We note that the buoyant effect may 
also drive the dewetting of air film. However when compared with the surface tension, 
the buoyant effect is completely negligible during the contraction period studied here: 
Ap • V • g/A ~ O.lPa <C a/5 ~ 10 3 Pa. Therefore we neglect the buoyant effect and use 
the surface tension as the solely driving force: 



Here v(t) is the instantaneous contraction velocity at the edge, and Ci, C2 are pre- 
factors of order one. This expression relates the contraction speed v(t) to the rim size 
5(t). From v(t), we can further calculate the pattern's radius, L(t) (see Fig. 5(a)), via the 
relationship v(t) = — dL(t)/dt. Thus we can obtain L(t) by solving for v(t) in Eq J3.3l and 
then integrate it: 



Therefore, L(t) can be calculated from the knowledge of S(r). We obtain the function 
6(r) by first measuring S experimentally and then fitting the data with high order poly- 
nomials, as shown in Fig. 5(b). On the other hand, L(t) can be independently measured 
from our images. To directly measure L(t), we make the best circular fit for the edge of 
the pattern, and obtain L(t) from the radius of the circle. The results of two approaches 
are compared in Fig. 5(c) for different liquids. The calculated curves match the measured 
symbols quite well, confirming the validity of the expression in Eq J3.31 In addition, the 
two fitting parameters, Ci = 0.63 and C2 = 1.54, are universally valid for all liquids, and 
are indeed of order unity. The air film contracts the fastest in H2O, due to the large 
surface tension and low viscosity. It contracts slower in oil- 1.04, because of the reduced 
surface tension. The motion is the slowest in oil-9.30, due to the effects of both low 
surface tension and high viscosity. In conclusion, the contraction is determined by the 
complex interaction between surface tension, inertia, and viscous effects, as described by 
Eq ]3.31 Because the dynamic contact angle during contraction is different from the static 
value and can not be measured directly, its role for air film contraction is still unclear. 



4. Conclusions 

In this work, we couple high-speed photography with optical interference to study the 
formation and evolution of air entrapment in liquid-solid impacts. We find a compressed 
air film formed before liquid touches solid, with an internal pressure significantly higher 
than the atmospheric value. After contact, the air film vertically expands at the edge, 
reducing the pressure within a few tens of microseconds and making a rim much thicker 
than the interior. This thick-rimmed structure subsequently contracts into an air bubble, 
driven by the complex interaction between surface tension, inertia, and viscous effects. 
Our investigation provides explicit information for the initial impact pressure as well as 
the detailed profile transformation during the subsequent evolution. This knowledge may 
benefit impact related applications such as surface coating and splash control. 



(3.3) 




(3.4) 
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